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ABSTRACT 

We show that by directly correlating the cluster kinetic Sunyaev Zeldovich (KSZ) flux, the cluster 
peculiar velocity power spectrum can be measured to ~ 10% accuracy by future large sky coverage 
KSZ surveys. This method is almost free of systemics entangled in the usual velocity inversion method. 
The direct correlation brings extra information of density and velocity clustering. We utilize these 
information to construct two indicators of the Hubble constant and comoving angular distance and 
propose a novel method to constrain cosmology. 

Subject headings: cosmology: large scale structure: theory-cosmic microwave background 



1. INTRODUCTION 

The large scale peculiar velocity field, arising from pri- 
mordial density perturbation, is a fair tracer of the large 
scale structure of the universe and is of great importance 
to constrain cosmology and the nature of gravity. 
But the huge Hubble flow contamination makes its 
measurement extremely difficult, thou gh several impor- 
tant p r ogresses have b e en ma de (e.g. ID avis fc Peebles! 



19831k iFeldman et"aT] [|2003ft h The cluster kinetic 
Sunyaev Zeldovich (KSZ) effect, which directly mea- 
sures the cluster peculiar momentum with respect to 
the CMB, opens a new wi ndow for the measurement 
of the peculiar velocity v„ llHaehnelt fc TegmarkHl996t 
Kashlmskv fc Atrio-Barandelal 120001: lAghlfflTmetHd] 
20011 lAtrio-Barandela et all 12004 iHolderl I2004D . But 



the direct inversion of v p requires extra measurements 
of the cluster Thomson optical depth r and cluster 
temperature. Contaminations in the thermal Sunyaev 
Zeldovich (TSZ) effect and KSZ measurements and 
inappropriate assumptions involved in this process put 
a systematic limit of ~ 200 km/ s to the inferred v n 
^IHH lAghanim et"aT] l200l iDiaferio et all 

A natural solution is to avoid the v p inversion, but 
use the KSZ flux, the direct observable. Contaminations 
to the cluster KSZ flux have different clustering proper- 
ties and can be applied to disentangle KSZ signal from 
contaminations. In the direct correlation estimator of 
the KSZ flux, many contaminations automatically vanish 
and most remaining contaminations can be subtracted. 
We will show that at z > 0.3, the systematics become 
sub-dominant and the statistical errors, at ~ 10% level 
for South Pole telescope (SPT 4 ), dominates. We fur- 
ther derive a novel method, which measures the Hubble 
constant and comoving distance and thus least model 
dependent, to constrain cosmology. Throughout this pa- 
per, we assume fi m = 0.3, tt\ = 1 — £7 m , <7s = 0-9 and 
adopt BBKS transfer function ijBardeen et al.lll98r^ . 



2. THE FLUX POWER SPECTRUM 

The KSZ cluster surveys directly measure the sum 
of cluster KSZ flux .Sksz and various contaminations, 
such as intracluster gas internal flow, radio and IR 
point sources, primary CMB, cosmic infrared background 
(CIB), etc. The signal is 

-Sksz = - ^ - ' T C mb(t)v p A£Ia = S W ov wo • (1) 

Here, (r) is the r averaged over the solid angle Af2^, 
5*100 is the KSZ flux assuming v p — lOOkm/s and 
viqo = Up/100km/s. Since at v ~ 217 Ghz, the non 
relativistic TSZ effect, which is one of the major con- 
taminations of the KSZ, vanishes, throughout this pa- 
per, we focus on this frequency. At v ~ 217Ghz, 
OB/dT = 540 Jy sr" 1 uK _1 . 

Optical follow up of KSZ surveys such as dark energy 
survey will measure cluster redshift z with uncertainty 
< 0.005. 5 The z information allows the measurement of 
3D correlation 



is(r) = {(Si-TA!U){Sj 
= {Sksz, % Sksz, j) + 



I An,)) 



(2) 



Here, I = yj Si/ Af2 j and Afi is the solid angle of each 
cluster. In the correlation estimator, many contamina- 
tions such as internal flow and instrumental noise vanish. 
The majority of remaining systematics can be subtracted 
directly and the residual systematics is generally (much) 
less than the statistical errors (32J- 

Since (SkszSksz} i s cluster number weighted and since 
uioo is correlated at large scale, the corresponding KSZ 
power spectrum Aj^ gz has contributions from both the 
velocity power spectrum A^ ioo and the power spectrum 

Ag vSv of (S1V1S2V2) , where 5 is the matter overdensity. 
We then have 



Al^(k) = ([ S 100 (m)^dm) Al iQQ {k) (3) 
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The photo-z of each galaxies has dispersion ~ 0.05. Clusters 
have > 100 galaxies and thus the determined z dispersion is < 
0.005." 
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Here, dn/dm is the cluster mass function. Throughout 
this paper, we assume a unity velocity bias and cluster 
number density bias (b n ) — 3. A 2 ioo is given by the linear 
theory prediction 



1 



,2 ^DltC 1 ! z ) 



(4) 



Here, k is in unit of /i/Mpc and (3 = (a/D)dD/da where 
13 is the linear density growth factor. E(a) is the evo- 
lution of Hubble constant normalized to E(a — 1) = 1. 
A^ M is the dark matter power spectrum (variance). We 
have assumed that two lines of sight are parallel to each 
other. For less than 10° angular separation, the accuracy 
of this assumption is better than ~ 1%. Because of the 
large bias (b n ), the (S1V1S2V2) term becomes dominant 
even in the linear regime at k > 0.06/i/Mpc. Because 
the nonlinear scale at z ~ 1 is k ~ 0.5/i/Mpc, we can 
still using the perturbation theory to predict 



x <5iwi<5 2 "2 



(k) = W(a)a 2 ^X 
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d s k 2 



-> -Ao M (fc, z)ct^(z) when fc > 0.15/i/Mpc 

Here, <x 2 100 is the one dimensional velocity dispersion 
in unit of (100km/s) 2 . v p is determined by large scale 
gravitational potential, so it decouples from small scale 
density fluctuation. Thus the last expression of Eq. 
holds even in the nonlinear regime, if substituting the 
two quantities with the corresponding nonlinear ones. 

The KSZ signal is amplified in the correlation, with 
respect to many contaminations due to three reasons. 
(1) Because of the k 2 denominator, A 2 , peaks at k ~ 
0.05/i/Mpc, while Ap M , which many contaminations fol- 
low, keeps decreasing toward large scales. (2) The extra 
^■S 1 v 1 s 2 v 2 term causes A 2 CSZ to keep increasing toward 
large k. (3) A 2 {SZ has much weaker z dependence at 
k < 0.03/i/Mpc (Fig. [TJ, comparing to A^ M . 

We have not considered the effect of redshift distortion 
in the above calculation. One can show that the redshift 
distortion is negligible at k < 0.1/i/Mpc in the A 2 part. 
Its e ffect to A 2 „,r can be dealt with by usual methods 
(e.g. iKaiserl (^987)). While the overall effect of the red- 
shift distortion is to increase the signal at k < 0.4/i/Mpc, 
it reduces the signal at k > 0.4/i/Mpc and thus makes 
the measurement of A 2 CSZ more difficult. Since we are 
mainly interested in the linear regime, for simplicity, we 
disregard this effect. 



3. NOISE ESTIMATION 

We target at SPT to estimate the accuracy of A|; gz 
measurement. We assume that SPT will cover the 217 
Ghz range and 4400 deg 2 with arc-minute resolution. We 
discuss the three dominant error sources, diffuse fore- 
grounds and backgrounds f i)3.1|l . sources associated with 
clusters f M3.2[l and cosmic variance and shot noise f i)3.3|) . 
We assume redshift bin size Az = 0.2 and k bin size 
Afc = 0.5k. 

3.1. Diffuse foregrounds and backgrounds 

The KSZ signal has typical AT KS z ^ 20(iK(t) /0.01. 
Primary CMB, which is indistinguishable from the KSZ 
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Fig. 1. — The normalized KSZ cluster flux power spectra at 
z = (black line) and z = 1 (blue) respectively. The dash lines 
are the contributions from the vv term and the dot lines are the 
contributions from the 5vSv term. The solid line is the sum of 
them. 



signal in frequency space, has intrinsic temperature fluc- 
tuation AT - lOOfjK. The CIB has mean . T ~ 20^K, 
if scaled to 217 Ghz ijFixsen et al.l 11998^ . A cluster 
KSZ filter can be applied to filter away the mean back- 
grounds while keeping the KSZ signal. This filter must 
strongly match the cluster KSZ profile while having zero 
integrated area. Since the angular size of clusters are 
generally several arc-minute, such filter naturally peaks 
at multipole / around several thousands and thus nat- 
urally filters away the dominant CMB signal, which 
concentrates at I < 1500. For such filters, at ~ 217 
Ghz, the galactic synchrotron, free-free foregrounds, the 
galactic dust emission, the radio background, the TSZ 
background are all negligi ble due to the i r frequency or 
scale d ependence (see, e.g. iWrightl lll998j) : lBennett et alJ 
(2003)). So we only discuss the contaminations of the 
primary CMB, CIB and background KSZ. 

The optimal filter can be constructed from the intr- 
acluster gas profile inferred from the TSZ survey. For 
simplicity, we choose the electron density profile as 
n e (r) oc (1 + r 2 /r 2 )^ 1 and a compensate Gaussian fil- 
ter W(l) = 6{l/l f ) 2 exp{-(l/l f ) 2 ). For these partic- 
ular choices, filtered KSZ temperature, ATksz, peaks 
at If ~ 1.1/0 C ~ 3800(1 /0 C ), where 8 C is the angular 
core radius, and the peak value is ATksz — ATksz — 
9wioo [(T)/0.01]fiK. For simplicity, we adopt this // and 
8 C = 0.4Mpc/h/x(z) to estimate the noise. Here, x( z ) is 
the comoving angular distance. 

The correlations of filtered backgrounds (with zero 
mean flux), originating from both their intrinsic 
correlations and cluster clustering, are ^fc(r) ~ 
[/ miow AQ A dn] 2 [l + (b n ) 2 ^ BM (r)}w b (9 - r/ X ), where 
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Fig. 2. — The dominant errors of the cluster KSZ power spectrum 
measurement. The solid lines are the errors caused by the cosmic 
variance and shot noise. The short dash lines are the systematic 
errors caused by the CIB and KSZ background. The long dash 
lines are the statistical errors caused by the primary CMB. The 
dot lines (which are visible only in the upper left panel) are the 
statistical errors caused by sources associated clusters. Systematic 
errors are sub-dominant at k < 1/i/Mpc. 



Wb(6) is the corresponding (filtered) background angu- 
lar correlation function. 

The CMB contamination is sub-tractable and does 
not cause systematic error. The same cluster survey 
measures both / Att A dn and (6„} 2 £dm, while the CMB 
Ci has been precisely measured b y WMAP and can be 
precisely predicted by CMBFAST (jSeliak fc Zaldarriagal 
1996). Thus the CMB contamination £cmb can be 
straightforwardly predicted and subtracted from the cor- 
relation estimator. Thus only the statistical errors 
caused by the CMB intrinsic fluctuation over cluster re- 
gions remains. 

The CIB and KSZ contaminations are in principle 
sub-tractable, too. But since both the amplitude and 
shape of CIB and KSZ power spectra are highly un- 
certain, we do not attempt to subtract their contribu- 
tion from the correlation estimator. The CIB power 
spectrum is C p m l 2 /(2ir) ~ (4^K) 2 (7/10 3 ) ' 7 (see, e.g. 
Zhang (2003)) and the kinetic SZ power spectrum is 

-,back-KSZ;2 l(0~A ~ (2.1 yK) 2 



iZhang et aW 2004). The 
upper limit of the fractional systematic error they cause 
is 



;6=CIB,KSZ 



[9((r>/0.01)^P 

i + <M 2 



A 2 , 



(6) 



,(k) + (b n ) 2 Al ViS2V2 ' 
which is at most several percent at k < 1/i/Mpc (Fig. |2J) 



The map filter does not filter away the contaminations 
from sources associated with clusters. At 217Ghz, the 
non-relativistic TSZ vanishes. But the relativistic cor- 
rection of cluster TSZ effect shifts the cross over point 
slightly to higher frequency and thus in principle intro- 
duces a residual thermal SZ signal in ~ 217Ghz band. 
We assume an effective AXtsz ~ l^K, or, ~ 1% of 
the TSZ at Rayleigh- Jeans regime. This gives a flux 
S T es-SZ ~ 540 Jy sr" 1 . The flux of radio sources and IR 
sources a ssociated with cluster i s ~ 10 3 Jysr _1 at 217 Ghz 
(see, e.g. lAghanim et alJ l)2004|) ). By multi-frequency in- 
formation and resolved source subtraction, one is likely 
able to subtract much of these contaminations. Rather 
conservatively, we assume that, at v ~ 217Ghz, the total 
flux contributed by those sources associated with clus- 
ters is less than ~ 5 x 10 3 Jysr^ 1 at z = 0.5 and scale it 
to other redshifts assuming no intrinsic luminosity evo- 
lution. 

The mean flux of these sources is subtracted in our 
estimator (Eq. [21 . Since the cluster thermal energy, IR 
and radio flux should be mainly determined by local pro- 
cesses, one can omit the possible large scale correlation 
of these quantities. Thus these sources do not cause sys- 
tematic errors. But since (5S 2 lst ) ^ 0, they do contribute 
to statistical errors. The fractional error they cause is 

-1/2 



(SS 2 lst ) (b n ) 2 A 2 BM 27r[Vk 2 AkY 
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(7) 



Here, V is the survey volume in the adopted redshift bin. 
It is reasonable to assume that (SS 2 ^) 1 ^ 2 < (Scist)- Thus 
the error caused by the sources associated with clusters 
is negligible at almost all scales and redshifts (Fig. [2J 

3.3. Cosmic variance and shot noises 

The signal intrinsic cosmic variance dominates at large 
scales. The amount of cluster is very limited, so the shot 
noise is large, even in the linear scales. The signal, the 
instrumental noise, the rms flux fluctuations of any con- 
taminations projected onto or associated with clusters, 
such as primary CMB, IR and radio sources, all con- 
tribute to the shot noise, whose power spectrum is 
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(8) 

Here, n is the mean number density of observed clusters, 
which can be calculated given the halo mass function, 
the survey specification and the gas model. For sim- 
plicity, we assume n(z) = 3 x 10 _5 /(1 + z) 3 (/i/Mpc) 3 . 
We only consider the listed 4 dominant sources, where 
ocmb ~ ctcib ~ ctksz ~ 20 and <7 clst ~ 5fiK. The 
beamed and filtered SPT instrumental noise has rms 
~ 1/j.K and is thus negligible. For SPT, the error caused 
by the cosmic variance and shot noise (~ 10%) dominates 
over all other errors. The systematic errors are (almost) 
always sub-dominant. In this sense, our method is opti- 
mal to measure the cluster peculiar velocity. For a future 
all sky survey, total error can be reduced to several per- 
cent level. 



3.2. Sources associated with clusters 



4. CONSTRAINING COSMOLOGY 
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Fig. 3. — The time variation of Fi(k = 0.015/i/Mpc (solid 
lines) and F 2 (k = 0.2h/Mpc, k a = 0.015/i/Mpc (dash lines) 
caused by wrong assumed cosmological parameters. Our fidu- 
cial model (black straight line) has f2 m = 0.3 and Q\ = 0.7. 
Blue, magenta, red and black curves correspond to (f! m ,f2A) 
; I M 0.0), (0.3, 0.0), (0.5, 0.5) and (O.f , 0.9), respectively. 



There are several ways to constrain cosmology from 
the KSZ observations, such as directly comparing the 
theoretical prediction of Eq. [3] with observations. Ex- 
tra information contained in the correlation allows us 
to develop a least model dependent method. Defining 
r)(k) = A^ sz (fc)/A^(fc), utilizing the asymptotic be- 
havior that n(k -+ 0) -» {S W0 ) 2 f3 2 E 2 a 2 /{b n ) 2 3k 2 and 
r)(k — » oo) — > {Sioo) 2 a 2 , we obtain two observables, 
Fi(k) = A|- sz (fc)/?7(fc — > oo) = Al(k)/a 2 in the regime 
where k < 0.03/i/Mpc, and F 2 (k,k ) = Al(k)/[g(k -> 

00)/rj(k o -> 0)] = A2 )M (fc)/[3CT2fc2/ /3 2 £ ;2 a 2] ^ thfl linear 

regime where fc < 0.3/i/Mpc. Fi^ should be z inde- 
pendent. Redshift distortion and possible velocity bias 
do not change this characteristic behavior. But if the 
cosmology one choose is different from the real cosmol- 
ogy, k w /k r , the ratio of determined wave- vector k w with 
respect to the real k r , will evolve with z. Because of 



the k dependence of Fi{k) and F 2 (k), a wrong cosmol- 
ogy will cause the determined F± and F 2 to vary with 
z. This behavior can be applied to constrain cosmology. 
One can show that kn w /k\\ r = E w {z)/E r {z) = gu(z) 
and kj_ >w /kx tr = Xr(z)/Xw(z) = 9±{z). Any mea- 
sured quantity ^(k) are then the corresponding quan- 
tity averaged over all configurations, namely, (^(A;, it)) = 



J\ * k^u 2 /g 2 + (l-u 2 )/gl du/2. Wc 



obta 



F hw (z) , (A 2 v (k,u)) 



Fi >w (Zq) 

F2, w (z) , 



Al(k) 
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(A 2 DM (k,u)) 
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Fig. y] shows the z dependence oiFi 2 {z) for various as- 
sumed cosmologies. This variation can reach ~ 20% and 
thus should be detected by SPT. A future all sky survey 
would measure this variation to better than ~ 5% accu- 
racy and should provide a tight and independent cosmol- 
ogy check. The nonlinearity correction, the parallel line 
of sight approximation, etc. may cause non-negligible 
systematics. But these effects are straightforward to pre- 
dict and correct, so we postpone the discussion of these 
effects in this paper. 



5. CONCLUSION 

We present a new method to measure the peculiar ve- 
locity power spectrum A 2 , by directly correlating the KSZ 
flux. The cluster peculiar velocity signal is amplified in 
the direct correlation. Many systematics involved in the 
usual v p inversion disappear and the majority of remain- 
ing systematics can be easily subtracted. The correlation 
method is optimal to measure A 2 ,, in the sense that sta- 
tistical error, dominates over systematics at almost all 
k and z range. We estimate that SPT can measure A 2 
to ~ 10% accuracy and a future all sky survey can im- 
prove this measurement by a factor of several. The KSZ 
flux correlation contains extra information on the veloc- 
ity dispersion. This extra information helps to construct 
two independent observables as indicators of Hubble con- 
stant and comoving angular distance, with only minimal 
amount of assumptions. These observables should put 
independent constraints on cosmology. 
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